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Abstract 

We calculate the area of a marginally trapped surface formed by a head-on collision of 
gravitational shock waves in AdSz). We use this to obtain a lower bound on the entropy 
produced after the collision. A comparison to entropy production in heavy ion collisions is 
included. We also discuss an 0{D — 2) remnant of conformal symmetry which is present 
in a class of gravitational shock wave collisions in AdS/j and which might be approximately 
realized (with D = 5) in central heavy-ion collisions. 
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1 Introduction and summary 

Relativistic heavy ion collisions produce a lot of entropy. Consider for example gold ions col- 
liding with y/sNN = 200 GeV. A head-on collision (usually described as "central" ) produces 
about 5000 charged tracks: see for example [T]. The total entropy may be roughly estimated 
as 

5 ^ r.SA^ehargcd = 38000 . (1) 

In section we explain where the factor of 7.5 comes from. 

The main aim of this paper is to inquire how well one can understand total entropy 
production in a heavy ion collision in terms of a dual black hole description. Ideally, we 
would like to construct colliding nuclei in a holographic dual to QCD. When the duals of 
the nuclei collide in the bulk, a black hole should form, signifying the formation of a quark- 
gluon-plasma. While a holographic dual to QCD is unavailable, it was suggested early on 
[21, [3] that an analogy should exist between colliding heavy ions and colliding gravitational 
shock waves in anti-de Sitter space. Subsequent related work on collisions in AdSs includes 
[H O E]. In the next few paragraphs, we will summarize an entropy estimate based on 
colliding gravitational shock waves which gives a result surprisingly close to ([T]). 
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The line element for two identical head-on shock waves propagating toward one another 
in AdSs is 

ds^ = ^ \-dudv + idx^f + idx'Y + dz^] + x^ z) \5{u)du^ + 5{v)dv^] , (2) 

^2 L J ^ L J 

where we have introduced the light-cone coordinates 



u = t — X 



3 _ + I ^3 



t + x^ (3) 



and have assumed that m < or f < 0. A simple shock-wave geometry in anti-de Sitter 
space can be obtained by boosting a black hole solution. As we will explain in section [3l for 
such a shock wave, the function z^, z) in ([2]) is given by 



^ ' L ,/^^iTTq) 



(5) 



where 

(xy + (x^)^ + (z-Ly 

4il 

and E is the total energy of the shock wave. In the rest of this introductory discussion we 
focus on shock waves of the form (jlj). Extensions to more general shock waves in various 
dimensions can be found in sections 13.21 and 13.31 

The metric ([2]) has singularities at m = g = and v = q = where Einstein's equations 
apply only in a distributional sense. These singularities merely signal the presence of pointlike 
massless particles of energy E, remnants of the boosted black hole. These singularities could 
be smoothed out by replacing each massless particle by a continuous cloud of massless 
particles with the same total energy. In [7], point-like sources for shocks propagating in a 
fiat-space background were replaced by wave-packets. The geometry ([2]) describes a head-on 
collision because the massless particles are located at the same position in the transverse 
space parameterized by x^, x"^, and z. 

The reason we must assume u < or < in ([2]) is that the two shocks collide at 
u = V = 0, and in the future light- volume of that event, little is known about the geometry 
(see however [6].) Assuming a black hole is formed after the collision, there is a standard 
method [SI [HI [IHl [H] for computing a lower bound on the entropy S of the black hole: 

Q ^ Q ^trapped /„n 

•J ^ -Jtrapped = , (.Oj 
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where ^trapped is the area of the trapped surface: that is, a surface whose null normals all 
propagate inward. The inequality (jSj) is based on the expectation that trapped surfaces must 
lie behind an event horizon. To our knowledge, ([6]) has not been rigorously demonstrated in 
anti-de Sitter space. It is related to singularity theorems [12], cosmic censorship (for a review 
see [13]) and the area theorem, which is usually proven on the assumption that spacetime 
is asymptotically fiat; see however [TH [T51 UHl [IZ]. Instead of attempting to clarify the 
conditions under which ([6]) must hold, we will make the working assumption that it does 
hold for the collisions we discuss. 

When they exist, trapped surfaces are highly non- unique. But in the case of head- 
on collisions in flat space there is a standard choice of such surfaces [HI El [101 E] which 
are easily obtained due to the symmetries of the configuration: head-on collisions preserve 
rotational symmetry around the axis of motion of the massless particles, 0(2) in ci = 4. 
The standard trapped surface preserves this symmetry too. In the case we're considering, 
the metric ([2D possesses an 0(3) symmetry which acts on x^, x"^, and z but preserves q. 
It is a remnant of the 0(4,2) symmetry of AdSs. We explain this symmetry more fully in 
section [3l and in section [H we construct an 0(3)-symmetric trapped surface (more precisely, 
a marginally trapped surface) which is an obvious adaptation of the standard one in fiat 
space. The marginally trapped surface we find comprises two halves. Si and S2, which are 
matched along a co-dimension three "curve" C. This is depicted in figure [H C lies in a three- 
dimensional slice of AdSs whose internal geometry is the hyperbolic space H3. Because 
of the 0(3) symmetry, C must be a two-sphere located at some constant value qc of the 
0(3)-invariant variable q. 

As we will see in section [H the area of the trapped surface depends on the energy E of 
the configuration, and one can obtain a relation between S'trappcd and E through qc- When 
gc ^ 1, this relation takes the form 

„ _ 4L2 3 _ AttL' 2 

~ ~P^1C "Jtrapped ~ 1c ' ^ ) 

from which we can immediately extract 

/ r3\ 1/3 

Strapped ^ ^ {2ELf/' . (8) 



G 



To obtain a numerical value for Strapped, we must evidently select values for the dimensionless 
quantities L^/G^ and EL. To choose L'^/G^, consider the translationally invariant AdSg- 
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Figure 1: A projection of the marginally trapped surface that we use onto a fixed time slice 
of the AdSs geometry. The size of the trapped surface is controlled by the energy of the 
massless particles that generate the shock waves. These particles are shown as dark blue 
dots. 



Schwarzschild solution: 



z 



According to [18], the energy density is 



16 Gs 

On the other hand, lattice calculationaH show that 



l-^\dt' + dx^ , 



(9) 



(10) 



for l.2Tc <T <2Tc 



and that rises slowly above this range. We choose 

16 

X 11 ^ 1.9 



G. 37r3 



;ii) 



(12) 



in order to make the black hole equation of state ffTOl) match ffTTl) . Since we have not specified 
a compact manifold, we need not assume that the AdSs background is dual to SU{N) Af = 4 



^We took the value quoted in (fTTjl from Figure 1 of [19]. See e.g. [20] for a more comprehensive account. 
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super- Yang-Mills. If we did, (fT2!) would imply that N ^ 2. Instead, we are assuming that 
the background is an approximate dual to real-world QCD above the confinement transition, 
or to a theory which is sufficiently close to real-world QCD to make numerical comparisons 
meaningful. Alternatively, we are assuming that the dual of the AdSs background captures 
enough features of real world QCD (above the confinement transition) to make this numerical 
comparison meaningful. In any case, loop effects on the gravity side are suppressed only by 
powers of G^/L^, so according to (fT2|) they are not very suppressed. Also, a' corrections could 
be signiffcant. Thus, all our calculations are to be understood as leading-order estimates. 

To choose a reasonable value of EL, we have to know a little more about the holographic 
dual of a shock wave. As we explain in section [31 the expectation value of the gauge theory 
stress tensor for the right-moving shock is 

1 IL'^E 
(T^Jx)) = lim —<l>(x\x^, z)6(u) = ^6(u) , (13) 

\ uuK )l 4^C'5^-0^3 V ^ } \ } (L2 + (xl)2 + (x2)2)3 ^ ^' ^ ^ 

with all other components vanishing when one uses the coordinate system (n, t', x^, x^). 
Evidently, E is the total energy in the gauge theory. For gold-gold collisions, ^/sNN = 
200 GeV means E = -Eboam = 19.7 TeV. L is the rms radius of the transverse energy 
distribution in (fT3l) . Because we're comparing the dual of the shock wave to a boosted gold 
nucleus, an obvious approach is to set L equal to the rms transverse radius of the nucleons. 
Using a Woods-Saxon proffle for the nuclear density (see for example [211 122]), one obtains 
an rms transverse radius L ^ 4.3 fm. So we estimate 

EL^4.3xlO^ (14) 

Putting dHD, dnD, and dUD together, we ffnd 

S > ^trapped ~ 35000 ( ] . (15) 

_ trapped V200GeVy 

In figure [2] we have plotted the dependence of the entropy bound (fTS!) on the energy, 
together with the data from PHOBOS [23]. It is encouraging that the estimate (fT5|) for 
'S'trappod is just 10% bclow the phenomenological estimate ([T|) at ^/sNN = 200 GeV. According 
to [23], once we use rj/s = l/47r, this is roughly the amount of entropy required to fit a 
thermalization time of 1 fm/c. However, the scaling Strapped oc s](^ implied by f|T5l) differs 
from the observed scaling, which is closer to the dependence S oc s^^f^. As observed in [26] . 
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Figure 2: A plot of the total number of charged particles vs. energy. The data points were 
taken from table II of the PHOBOS results [S^. We show in red the region consistent 
with the bound ffT^ obtained via the gauge-string duality, using point-sourced shocks and 
estimates described in the text, and assuming the bound The blue curve corresponds to 
the prediction of the Landau model 

the latter dependence, predicted by the Landau model [23] seems to hold over a strikingly 
large range of energies. Put differently, the inequality in fll5|] is consistent with all heavy-ion 
collision data to date, but for energies only slightly above RHIC energies, (I15p predicts a 
faster increase of entropy than is generally expected. 

At the LHC, y/s^N will be 5.5 TeV for lead-lead collisions. Inserting this value into (fTSjl . 
and making minor corrections for the differences between lead and golcj^l one finds 



S > 5'trappcd corresponds to A'chargcd > 45000 if we continue to use ([T]). The lower bound 
on the entropy f|T6l) exceeds the prediction of the Landau model, S ~ 2.1 x 10^, by a 
factor of about 1.6. Calculations based on the Color Glass Condensate tend to predict 
lower multiplicities: for example, from figure 5 of |28], one may read off the prediction 
A'chargcd ~ 22000, about a factor of 2 below the estimate from ([TB|) : see also 
We see three main ways in which ( HM could fail: 

1. Using the gauge-string duality to describe entropy production may cause us to mis- 

^For an introduction to the Landau model, see section [2T3l for a review, see [27] . 

= 208 for lead, so i?bcam = 570 TeV; L = 4.4 fm from the rms radius of lead, resulting in EL « 1.3 x lO"^. 



'trapped 



,4 X lo^ 



(16) 
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represent perturbative aspects of the early stages of the colhsion. This is because our 
use of the gauge-string duahty rehes on the supergravity approximation, which is the 
leading order description of a strong coupling expansion as well as a expansion. 
Our methods appear to offer no access to perturbative physics. Perturbative QCD is 
expected to characterize the early stages of LHC collisions more cleanly than it does 
RHIC collisions, and it may be that our methods are correspondingly less applicable 
at LHC than at RHIC. 

2. As we will see in section 13. 4[ there is a whole family of AdSs shock waves with the 
same (Tuu), presumably distinguished by higher point functions of Tuu- The trapped 
surface depends on which of these shock waves we pick. It is easy to lower S'trappcd 
by spreading the shock wave out over the transverse H-^ in AdSs. So although (fT5|) at 
first looks highly predictive, and easily falsifiable at energies significantly higher than 
RHIC scales, it is in fact possible to accommodate slower growth of total entropy with 
beam energy. We discuss this further in section [51 

3. The bound ([6]) could fail, even for standard Einstein gravity in AdSs. 

The rest of this paper is organized as follows. In section [21 we review phenomenological 
estimates of the entropy produced in a heavy ion collision, with the aim of justifying and 
briefly summarizing the dependence on beam energy. In section [3l we review the construction 
of shock waves in AdS/j, with particular attention to the 0{D — 1) symmetry preserved by 
head-on collisions of the simplest shock wave constructions. These symmetries might be 
approximately realized in central heavy-ion collisions even if gauge-string methods fail to 
give a quantitatively accurate account of entropy production. Aside from discussing these 
symmetries, our purpose in sections [2 and [31 is mostly to gather together well-known facts 
from the literature. Our main calculations are in section [H where we compute the shape of 
marginally trapped surfaces. We end with a discussion in section [51 

2 Phenomenological estimates of the entropy 

In order to evaluate the entropy S produced in a heavy-ion collision, one needs a method to 
relate the entropy to a quantity which can be measured: the number of charged particles, 
A'chargod- In l2.1l and l2.2l wc review two such methods. The first uses the framework of Bjorken 
flow pO] . The other, described in section 12. 2[ relies on phase space estimates to evaluate the 
ratio ^'/A'^charged after hadronization. Both of these sections largely follow [31] . 
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While sections 12.11 and 12.21 allow an evaluation of the entropy via the measured number 
of charged particles, in section 12.31 we estimate the entropy from the size and shape of 
the colliding nuclei — or, more precisely, the size and shape of the parts of the nuclei that 
participate significantly in the collision, and the beam energy per nucleon. This last estimate 
is based on the Landau model [21]. 

Non-specialists may appreciate the reminder that s at a? is the Mandelstam variable for a 
pair of nucleons, one from each nucleus. When the beam energy is 100 GeV per nucleon, 
^/sNN = 200 GeV. Because gold has 197 nucleons, the total center of mass energy is 39.4 TeV. 
It is also good to know that the rapidity of a particle emerging from the collision is y = 
tanh.~^ Pz/ E, whereas pseudo-rapidity is 77 = tanh.~^ pz/p = tanh~^cos6', where 6 is the 
angle from the beamline. 



2.1 Entropy estimate from Bjorken flow 

In this section (as well as in parts of section 12. 2p we estimate the total entropy produced in 
the collision by assuming that the entropy per charged particle changes only slightly with 
rapidity: thus 

S dS/dy 



(17) 

mid— rapidity 



-^charged '^-'^chargcd/'^y 

Both A'chargcd and dNchaxgcdl dy at mid-rapidity are directly measured, so we only need to 
estimate dS/dy at mid-rapidity. To do this, we follow [31] and consider Bjorken's boost- 
invariant treatment of a collision. 

One of the main relations emerging from Bjorken's treatment is 



^form^ f (^form) 



dy 



where e is the energy density, A is the cross-sectional area of the participating nucleons, rform 
is the formation time, and Et is the transverse energy of a particle, defined as E sin 6 where 
E is the total energy and 6 is the angle from the beamlinejj The entropy may be expected 
to follow a similar relation: 

dS 

Tform^ S(rform) = -y- ■ (19) 

dy 

If the QGP is a thermalized plasma at a time Tform (which may not be true but provides a 



''Sometimes the definition of Et is varied across particle species by adding some multiple of the rest mass: 
see for example [32] . 
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rough estimate), then assuming conformal invariance one has 

s = --. (20) 

By solving (fT8|) for e(rform) and (fT9|) for s(rform)) and then plugging the resulting expressions 
into ( l20l) . one arrives at 



dy 3T dy ' ^ ' 

The quantities in (pTj) are all to be evaluated at Tfornu but for a rough estimate one may use 

dy drj 

for central gold-gold collisions at ^/sNN = 200 GeV [I9]lf| To find the temperature, one may 
resort again to f|T8l) together with the approximation ffTTl) of the lattice equation of state. 
Choosing the nominal values Tform = Ifm and A = 120 fm^ leads to T = 240 MeV. Then 
(ETD together with ([22]) give 

c?g(rform) ^ 3300^ (23) 
dy 

The number of charged tracks per unit rapidity for a central collision is 

^ 660 (24) 

dy 

near mid-rapidity (see for example P). Once again as a rough estimate, let's use 

dS(fina\) dS^Tform) 



dy dy 
Putting (!23|) - (!25|) together, we arrive at 

C^»S'(Tform) ^ -^charged 

dy dy 

Using ( ITTl) leads to 



(25) 



(26) 



S ^ SiVchargcd ~ 25000 , (27) 
where we recalled that A'chargcd ~ 5000 in a central collision. It's important to bear in mind 



^V*NN = 200 GeV means that a pair of nuclcons, one from each nucleus, together have Mandelstam 
s = (200 GeV)^. That means the beam energy is 100 GeV per nucleon, and because gold has 197 nucleons, 
the total center of mass energy is 39 TeV. 
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that many approximations were used in arriving at (fT7|) . so it should be regarded only as a 
first attempt. Some refinements were outlined in |31j . 



2.2 Entropy estimates from phase space density 

Phase space estimates of entropy start from the expressions 

(Px d?p 



(27r)= 

where /j = fi{x,p) is the phase space density for each spin polarization of each hadronic 
species, and Sj = 1 for bosons and —1 for fermions. The number of hadrons is roughly 
l^chargcd- Usiug fl2Sl) mcaus that one is ignoring interactions among hadrons, so it should 
apply, in some approximation, after hadronization. One line of thought [531 Ell ES] is to use 
the equilibrium expressions 

1 

fi = r-, — , (29) 



set T = 170 MeV, and run the sums in (l28l) over all established hadron resonances. The 
result is 

5/iV = 5.15. (30) 
Applying fl30l) to the total entropy in a heavy ion collision gives 

5/iVcharged = 7.7. (31) 

The estimate fl3Tl) is at best approximate, because chemical potentials for quarks become 
significant at forward rapidities. Replacing S and Ai'chargcd by dS/dy and dNchargcd/ dy in 
f l5T]) would improve the reliability of the estimate. But fl5T]) is also approximate because the 
system is not really an equilibrated gas of nearly free hadrons at T = 170 MeV; rather, it is 
at roughly this temperature that the quark-gluon plasma hadronizes. 

A more data-driven approach was taken in [HT]: instead of assuming (1291) . experimental 
results for single-particle yields and two-particle interferometry were used to estimate the fi. 
For central collisions, and at mid-rapidity, one finds from this approach the result 

dS 

— = 4451 at = 130 GeV . (32) 

dy 
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Combining (132|) with 



(see for example [23]) gives the ratic|^ 



^^/^^ 7.2. (34) 



We arrived at the figure ^ simply by using 5'/A''charged = 7.5, an average of (|3T|) and ([3 



2.3 Entropy estimates from immediate equilibration 

The Landau model of particle production in high-energy collisions [23] assumes that hydro- 
dynamics is valid starting from the moment that the colliding nuclei completely overlap. It 
also assumes a conformal equation of state, p = e/3. The validity of hydrodynamics depends 
on local thermodynamic equilibrium and a mean free path that is short compared to the ex- 
tent of the medium. Total overlap occurs about 0.13 fm after the nuclei first start to interact 
in a central gold-gold collision at ^/SNN = 200 GeV. It doesn't seem reasonable to assume 
that hydro is valid at such an early time, so it is surprising how well the model works in 
describing aspects of the bulk fiow, in particular the particle distribution in rapidity. 

The entropy is easy to estimate at Toveriap = 0.13 fm. We should ignore the nucleons 
which do not interact: this includes a good fraction of the ones in the outer skin, or corona, 
of the gold nucleus. In a central collision (more precisely, in the 5% of collisions that are the 
most central) a typical number of participating nucleons is A^part = 350, so the total energy 
of participating nuclei is 

^^^^ = ^Vl|^ = 35TeV. (35) 

In the rest frame of one nucleus, its participants occupy a roughly spherical region of radius 
6.5 fm, which we will assume to have uniform density. Let's denote the volume of this sphere 
by V. In the lab frame, this sphere is Lorentz fiattened by a factor 7 = ^/snn /2mp, where 
mp = 0.938 GeV is the mass of a proton. Thus the energy density at the moment of overlap 

^Thc STAR collaboration has published a result corresponding to dNchmgcd /dy ~ 580, which would result 
in a value 7.7 in (j34p . An even higher value, 8.5, can be read off from estimates in [34]; however there seems 
to be some possible confusion about dNchargcd/dy versus dNch_a,rgcd/dr]. We have used a common mid-rapidity 
conversion factor, dNchargcd/dy ~ l.ldNchargcd/ dr/, to pass from results quoted in terms of pseudo-rapidity 
densities to rapidity densities. 
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IS 



lEtot 



3300 GeV/W . 



(36) 



Using (ITT]) , the corresponding temperature is 



T 



1/4 



1200 MeV , 



(37) 



and the entropy is 



S 



'tot 



^(4/,m,\/iV3^rtS^^)i/^ ^ 38000, 



(38) 



3 T 



which is fortuitously close to ([T]). If one started instead by assuming that all the nucleons 
participate and that the radius is 7fm, the entropy estimate would increase to 44000. The 
usual assumption in the Landau model is that subsequent expansion is isentropic. 

To arrive at the figure S ~ 2.1 x 10^ for entropy production at the LHC, quoted below 
f fT6|) . we used fl38l) with A^part = 368 (scaled up from the number for gold in linear proportion 
to the atomic number), R = 6.6 fm (scaled up from the number for gold in proportion to the 
cube root of the atomic number), and the same value = 11 as beforej^ 

3 Shock waves in anti-de Sitter space 

Gravitational shock waves are well studied, both in R"^^^'^ and AdSz): see for example 
[36l[371[38l[39lll0lllllll3ll3lll2ll5]. The simplest of them can be constructed in two 
equivalent ways. One is to boost a black hole in AdS/j to a velocity approaching the speed 
of light, while at the same time decreasing the mass of the black hole in such a way that the 
energy remains fixed. We describe this construction in section [3?T] for the special case D = 5. 
Alternatively, one can start off with a pointlike, massless particle traveling in AdS d and show 
that it back-reacts on the metric in such a way as to produce a shock-wave discontinuity. We 
describe this construction for arbitrary D in section 13. 2[ Other types of shocks obtained by 
sourcing the metric with appropriate matter are given in section 13. 3[ Our goal is to relate 
collisions of shock waves to collisions of heavy ions, and to this end it is useful to have the 
energy density dual to the colliding shocks. We compute this in section 13. 4[ 

''A fractionally higher value, say /* — 12, might be closer to lattice values, but it doesn't make a difference 
at the level of accuracy we have quoted. 
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3.1 Constructing the simplest shock wave geometry 

Our starting point is the global AdSs-Schwarzschild (GAdSBH) metric, 

ds' = -fdr' + ^ + p'dnl / = l + ^-4' (39) 

where the parameter po can be related to the ADM mass of the black hole by 

M=^/o- (40) 

Since we are working in global coordinates, the boundary theory has topology x R. 
Working in a coordinate system which covers only the Poincare wedge of AdS corresponds 
to putting the boundary theory on R^'^. 

To make this more precise, recall that AdSs is the universal cover of the five-dimensional 
hyperboloid 

-{X-'f - {Xy + {X'Y + {Xy + + (X^)2 = (41) 

in R^'2. The metric of AdSs, which is given by (139|) with po = 0, is also the metric induced 
on the hyperboloid from the standard flat metric on R^'^. The X'^'^ coordinates are related 
to the global coordinates (r, p, Q^) in fl39|l as follows: 

X-^ = + L2 cos y X° = Vp2 + L2 siu y X' = pQ' , (42) 

where is a unit vector in R^, which is to say a point on . But r runs from — oo to oo in 
AdSs, whereas r = and r = 2'kL are identified on the hyperboloid. Thus, the coordinates 
X*^ are more fit to describe the hyperboloid of which AdSs is the covering space, while in 
the (r, p, f2*) coordinate system, AdSs can be thought of as a cylinder with boundary S''^ x R 
if we conformally compactify in the p direction. A more detailed discussion can be found, for 
example, in [16]. The Poincare coordinates (t, x^, x^, x^, z) are related to the X^^ coordinates 
by 

2^ ' > ' (43) 

X^ = L- X'='-(-l + - +^ 



~2 



z 2 

The actual metric obtained when transforming fl39l) to the Poincare patch is somewhat 
complicated, and we shall not write it explicitly here. 
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Figure 3: The hyperboloid whose covering space is AdSs, with the transverse coordinates X^, 
X^, and X^ suppressed. The closed green curve is the trajectory of a massive test particle. 
When the particle is infinitely boosted, so that X^ = X^, the trajectory degenerates into 
the two blue lines. 

We wish to boost this black hole to the speed of light while decreasing its mass and 
keeping its energy constant. The method we use is similar to the one explained in |38j . 
Following [36] we expect that this boost will give us the gravitational field around a massless 
particle moving in AdSs. Thus, we choose a boost that will take a stationary particle to one 
moving at a highly relativistic speed. A massive test particle follows a closed trajectory which 
is described, in terms of the coordinates, by the intersection of the hyperboloid (HTj) with 
the plane X^ = (3X^ with fixed X\ X^, and X^. A convenient choice is X^ = X^ = X'^ = 0. 
In figure [3] we have shown one such trajectory. Once we take the mass of the particle to zero, 
the trajectory of the test particle degenerates into two straight lines (straight both in the 
sense of being geodesies on the hyperboloid and in the sense of the fiat metric of R^'^). This 
is also depicted in figure [31 In global coordinates, a massive particle starting at p = with 
some non-zero velocity returns to p = with the opposite velocity after a time r = irL, then 
continues to oscillate through the AdSs cylinder with period r = 27rL. When taking the 
lightlike limit of the trajectory, the oscillating motion of the massive particle deforms into a 
bouncing motion, going from one boundary to the other. Each leg takes a finite global time 
At = irL and corresponds to one line on the hyperboloid. 

To carry out the boost explicitly we note that the isometry of the hyperboloid fj^Tl) is the 
0(4,2) group of linear transformations preserving the quadratic form on the left hand side 
of (HTj) . The boost that we'll be interested in is an element of the 5*0 (1, 1) subgroup which 
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preserves X ^, X^, X^, X'^, and the quadratic form — (X°)^ + (X^)^: 

-/?X0 + X3 



X° = ^° X^ ^ X3 



(44) 



Writing M = E^Jl — fS"^ and taking /3 — >■ 1 with ii^ held fixed, the boosted GAdSBH metric 
fl39l) becomes 





\xy - (x-i)2" 


+ 


"(X0)2 + (X-l)2" 


(X-l)2 




(X0)2 + (X-l)2 


-L2 + (X0)2+(X-1)2 





(t/X° - rfX3^2 



Using 



f{x)dx 



which holds for any integrable function /, the /3 — >^ 1 limit of (l45l) becomes 



(45) 
(46) 



Wc2 _ ^.2 _ 

as — as^^g^ 



L2 - 2(X-i)2 + 2X-^y/{X-^y - 



L2^(X-i)2 -L2 
X e(X-^ - L)5(X° - X=^)(rfX° - rfX3)2 



(47) 



where 



e(x) = 



1 if X > 



(48) 



if .X < 0. 

In Poincare coordinates fH3|l . the line element P7|) reads 



(is^ = :^ r-dwdt; + (rfx^)2 + (dxY + ^^^'1 + -'^'(a;\ z)5{u)du^ , (49) 

where v, and $ are defined as in Q-Q. In checking the equivalence of (l47j) and (l49l) . it 
helps to note that 



X~'-L 
2L 



(50) 



A subtlety in the derivation above is the emergence of the factor of 6(X~^ ~ L). As 
we've explained earlier, the massless trajectory comprises of two straight lines describing a 
massless particle which goes back and forth from one boundary of the AdS cylinder to the 
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other. We choose to consider only one such leg, and this is what the factor Q{X~^ — L) does. 
Including the return journey of the particle would correspond to adding an additional term 
to p7|) identical to the one explicitly shown, but with X^^ — 

The metric (1491) has an 0(3) symmetry which is simpler to understand in the X^'^ co- 
ordinates: the boosted metric fHSj) doesn't depend on X^, X"^, or X^, except through the 
constraint fHTl) . which (after the boost) can be regarded as a way to determine X~^ in terms 
of X^, and {X^Y + (X^)^ + (X")^. The 0(3) symmetry is the one acting on the coor- 
dinates X^, X^, and X^ transverse to the particle's trajectory. To see the transverse space 
more clearly, we slice AdS5 at a definite value of X^ and impose X^ = X^. This gives a 
two-sheeted hyperboloid: 



-{X"'Y + {X'Y + {Xy + {Xy = -L^ . (51) 

The two disjoint sheets are related by X~^ —X^^. Each is a copy of the Euclidean 
hyperbolic space H^. The massless particle that we're interested in passes through the 
"center" of the upper sheet, at X~^ = L and X^ = X'^ = X'^ = 0. (We write "center" in 
quotes because is a homogeneous space.) The isometrics of H-^ form 0(3, 1)/Z2, and the 
0(3) of interest is the part of this group that preserves the point that the massless particle 
passes through. Note that the 0(3) symmetry we've found is not equivalent to rotational 
symmetry in the (x^, x^, x^) plane of the Poincare patch. Rather, because it acts non-trivially 
on X^, its generators include special conformal transformations. 



3.2 Shock wave metrics in AdS 



D 



The metric (jlSD can also be obtained by solving the Einstein equations in the presence of 
a lightlike particle. This alternative derivation is a little more efficient, and we will take 
advantage of this to generalize to D dimensions. Following [371 El 12] , one starts with an 
ansatz 

^2 / D-3 \ 

ds'^ = — { -dudv + y^(dx'Y + dz^ + (f)(x\ z)6(u)du^ . (52) 

^ V / 

If = 0, this is the metric of AdS/? in Poincare coordinates. Because is supposed to be 
the metric in the presence of matter, it should satisfy the appropriate Einstein equations: 

Rfiu — -j^QtiuR — ^~9iiu = SnGDJfiu , (53) 
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where J^i, is the bulk stress tensor, not to be confused with the boundary stress tensor 
For a massless particle with energy E, the only non-zero component of J^^, is 



D-3 



Juu = E 5{u)5{z -L)\[ 5{x') . (54) 

i=l 

It is straightforward to plug (152!) and (!54|) into the uu component of (!53|) and explicitly derive 

-^^^) ^ = -16vrGz5 E5{z -L)l[ 6ix') , (55) 



1=1 
where 

$ = -0 (56) 
z 

and 

i=l ^ ' 

is the laplacian on the Euclidean hyperbolic space whose line element is 



i=l 



Evidently, Hd-2 is the space transverse to the trajectory of a massless particle. If we 
introduced global coordinates X^'^ on AdSz), it would have a description entirely analogous 
to the one explained around (15T1) . For our present purposes, it is enough to introduce a 
subset of the global coordinates, as follows: 



L ■ 

= —x^ for i = 1 through — 3. 

z 



(59) 



ifD_2 is the upper sheet of the two-sheeted hyperboloid 

-{Vy + ^{Y'f = (60) 
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which has isometry group 0{D — 2, 1)/Z2E In analogy to (H2!) one may define 



yo ^ ^_ ^2 = for z = 1 through D - 2, (61) 

where 17* is a point on a unit S^~^. The metric (1581) can be re-expressed as 

dsl,-. = Y^^2+r'dnl_,. (62) 
We will be especially interested in the quantity 



i=l 



where the second equality can be checked using The last expression in (I63p shows that, 
up to an overall prefactor, q is the square of the chordal distance between the point {x\z) 
on H£,_2 and the special point {xl,z^) = {0,L) through which the massless particle passes. 
Chordal distance, by definition, is the distance in the embedding space R^~^'^ parameterized 
by the coordinates {Y^,Y^). In the coordinates introduced in (1611) . the special point is at 
r = 0, and 

-1 + a/1 + rVL2 ^- 

q = ^ — r = 2L^q{l + q). (64) 

Thus one may re-express 



L^ + 4,(l + ,).(7?,_3^ 



(65) 



The 0{D — 2)-symmetric solutions to (!55|) can be efficiently found by rewriting it in terms 
of q: 



where 



.^_3 _{D- 2)7r(^-2)/2 



Vol 5^-2 = ^ . (67) 

T{D/2) ^ ' 



^The isometry group of the full two-sheeted hypcrboloid is 0(D — 2, 1). The Z2 that one must divide out 
when considering a single sheet acts by sending Y'^ —Y'^. 
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The solution to (!66l) with the boundary condition that $(g) — as g — oo is 
It is easy to check that (jl]) is recovered by setting D = 



3.3 Other sources as shock waves 

Although gravitational shock waves are solutions of the full non-linear Einstein equations 
(in a distributional sense), two shocks moving in the same direction can be superposed: that 
is, if one starts with 

z 

as the first shock and 

ds^ = ds\as^ + —^2{x\ z)6{u - U2)du^ (70) 
as the second, then the superposed solution is 

ds'^ = ds\^SD + ~ [*i(^N ^)'^(^ - ^i) + ^2{x\ z)5{u - Us)] du^ . (71) 

If (l69l) and (170|) are sourced by massless point particles, then (!7T!l describes the back-reaction 
of the two massless particles together. More generally, we can consider a cloud of massless 
particles, all moving in the same direction, and then the metric is 

ds^ = ds\^So + -^(^*' «)^«^ • ( ''2) 

z 

The only non-trivial component of the Einstein equations is 

and it is straightforward to show that it takes the form 

Dho-2 -^^^) ^ = -iG^Gn^Juu . (74) 
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While it may be interesting to consider the case where F is non-zero over a range of u (see 
in this connection the recent work [IH]), let us restrict attention here to the case where 

F = z)S{u) and J„„ = -p{x\ z)5{u) . (75) 

z 

Then ( 1711) becomes 

D-2 



nHn^2 j^j * = -IQ^GdP ■ (76) 

Choosing p = E5{z—L) Y[f=i^ 5{x^) would return us to the point-sourced shock waves that we 
have focused on up until now, as can be seen from comparing fl76|) to fl55|) . The most general 
0{D — 2)-symmetric shock localized at n = corresponds to a source term p depending only 
on the chordal distance variable q defined in (ES]): then fITBl) becomes 

g(l + g)$" + ^(1 + HiD - 2)$' - {D - 2)$ = -IQtiGdL'p . (77) 

To solve (1771) . we follow the classic approach of first solving the homogenous equation and 
then using a Green's function to solve the general inhomogeneous equation. The solutions 
to the homogenous equation are 

$_(g) = l + 2g 

(78) 

<^+iq) = q'-''2FiiD - 2, D/2; D; -1/q) . 

Note that ^~{q) is the unique solution that remains finite at g = 0, and $+(q') is the unique 
solution that decays to zero at infinity. The Green's function G{q, go) satisfies 



1 



q{l + q)d' - -(1 + 2q){D - 2)d, -{D-2) 



2 



2 



Giq, go) = - [g(/7g^p-4)/2 ^(g - go) • (79) 



G(g, go) is uniquely specified by the requirement that when go > 0, G{q, go) should be finite 
at g = and should decay to as g — oo. Straightforward calculations lead to 

, SttGdL^ f'^'+(go)'^'-(g) forg<go 
G(g,go) = ^y^{ (80) 
^ [<l>_(go)$+(g) forg>go. 

The solution to the original problem (1771) is 

POO 

m = / dq, [go(l + go)](^-'^/' G(g, go)p(go) . (81) 
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Assuming that p has compact support, or else decays quickly enough at infinity, the asymp- 
totic behavior of $ near infinity is 

where we have defined 

/•oo 

E = 2^-3L^-2(Vo1 5^-3) / rfg + ^ 2g)p(g) 



Jhd-2 



(83) 



The power oi L/z in the second line is from the measure on Hd~2 associated with the metric 
fISS]) . The asymptotic expression fl52]) for $(5) coincides with the solution fl68p for a point- 
sourced shock wave. This amounts to a sort of shell theorem: an 0{D — 2)-symmetric cloud 
of massless particles gives rise to the same gravitational field, outside the cloud, as if the 
cloud were replaced by a single massless particle at its center with energy E. 



3.4 The gauge theory stress tensor of colhding shocks 

Before proceeding to calculate the trapped surface associated with the colliding shocks, we 
make an aside to discuss their dual boundary theory stress energy tensor. The holographic 
image of a shock wave geometry on four dimensional Minkowski space has a stress-energy 
tensor that can be found from the small z asymptotics of $. If 

ds^ = ds\^s, + ^^g^ndx'^dx^ , (84) 
where 5gmn ~ C>{z^) for small z and has no non-zero components with a z index, then |l9] 

^3 ^ 

{Tmn) = -r-pT 1™ -^Sgmn ■ (85) 

(This form of {Tmn) holds when the boundary metric is chosen to be —dt"^ + dx^. Conformal 
transformations on the boundary require a modification of ( l85i) .) We have used the notation 
Sgmn to represent the deviation of the metric from empty AdSs even though this deviation 
is not necessarily small. Applying the general rule (l85l) to the shock wave metric (H9|) . one 
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finds 

with all other components vanishing when one uses the coordinate system (m, t;, x^, x^). This 
same result may be obtained (as it essentially was in [H]) by first computing the stress tensor 
for the unboosted black hole and then applying the boost (144|) directly to {Tmn)- 

For D 7^ 5 the analysis is similar: the energy momentum tensor associated with the 
metric ( 1521) may be evaluated using ( l68l) and 



{D-l)L^+\ . 1 
= — -— hm -^dg^n , (87) 



the equivalent of ([HS]) [19]. We find that 



D-l 



TT 2 



where = Yl!i=i ^"i- ^^^o 

'"rf^-2a;(Too) =i?. (89) 



The profile ( IHTIl respects the 0(-D — 2) symmetry discussed at the end of section [321 

For the configurations discussed in I3.3[ we learn from comparing fl68|) to fl82|) that the 
dual expectation value for the gauge theory stress tensor must coincide with fl88|l . This 
illustrates a large ambiguity in the gravity representation of some given configuration of 
{Tmn)- Such an ambiguity should not be surprising since a state in the gauge theory is by no 
means completely specified by the one-point function of stress tensor. Knowledge of higher 
point functions of the stress tensor, and possibly of other gauge-invariant operators, would 
resolve such ambiguities. 

Going back to D = 5, we are eventually interested in collisions of nuclei. So we'd like 
to tune L and E in fl86l) to resemble the energy density of a boosted nucleus as closely as 
possible. For a gold nucleon at rest, the energy density of the nucleus can be read off of the 
Woods-Saxon number density n{x^ , , x^): 

e(x\x^x^) ^ mpn(x\x^x^) oc ^ ^ J^x\-R)/a ' ^'^^^ 
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where typical values for gold are [HI [22] 

R = 6.38 fm a = 0.535 fm . (91) 

To compare the energy density with (l86ll we need to boost (pOj) . Consider boosting stationary, 
pressureless dust: before the boost, 

(Too(^,a;^x^x3)) = e(x^x^x3), (92) 

with other components vanishing in the coordinate system {t,x^,x^,x^). After the boost, 
the non-zero components are 

\Too) {To3)\ 1 / 1 -1 



^T3o) {T,,)) l-P' \-l 1 

\/(3 -1 
-1 P 



e{x\x\{x''-pt)/^l^) 

(93) 

Pr{x^,x^)5{u) 



where 

PR{x\x^) = ^J==[dx\{x\x^,x^). (94) 

If the limit /? — > 1 is taken with p/j(x^,x^) held fixed (i.e. scaling e down by a factor of 
a/1 — /5^//5), then ( 193|) becomes simply 

(T„.)=Pii(x\x2)5(n). (95) 

Applying fl^ to fpl . we can find an expression for the uu component of the stress tensor 
of a nucleus. 

The resulting dependence of {Tuu) on the transverse position xt = \/ (x^)^ + (x^)^ is 
shown in figure HI together with the result (l86l) obtained from the simplest gravitational 
shock wave. The value of L for the gravitational shock wave was chosen so that the rms 
transverse radius (weighted by Pr{x^, x^)) matches to the same quantity computed using the 
boosted Woods-Saxon profile. Explicitly, for the Woods-Saxon profile, 

2, Jd^x [{x^f + {x^f]n{x\x\ x^) _2j d^x\x\y{l + e^^l-^)/'^) 



/rf3xn(xi,x2,x3) 3 /c/3a;/(l + e(l^l-«)/«) 

(96) 



Li5(-e«/'^) 



Li3(-e«/«) 
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Figure 4: Left: The dependence of (T„„) on transverse radius Xt, both for an infinitely 
boosted black hole in AdSs and for an infinitely boosted Woods-Saxon profile. (T„„) is 
proportional to S{u), and the quantities that we plot omit this singular factor. The normal- 
ization of the Woods-Saxon profile was chosen so that its maximum is 1. The normalization 
and width of (T„u) from the gravitational shock wave was chosen so that the integral and 
rms transverse radius match to the values extracted from the Woods-Saxon profile. Right: 
The area under the curves xt{Tuu) obtained from Woods-Saxon and AdS profiles are the 
same, indicating that the total energy is the same. 

where Li„ are polylogarithm functions. On the other hand, (x^) = for the AdS profile 
Pr{x^,x'^) oc [L^ + (x^)^ + (x^)^]""^. By plugging (19T|) into (p6l) and setting the result equal 
to L^, one obtains L ^ 4.3 fm, as quoted in section [TJ 

4 A marginally trapped surface for a head-on collision 

Once the two shocks collide we can no longer superpose the solutions for two single shocks. 
We assume that such a head-on collision will result in the creation of a black hole. A standard 
calculation in fiat space [H [H [TOl [11] is to estimate the area of the resulting black hole by 
constructing a particular trapped surface that lies on the t < parts of the u = and v = 
hyper surf aces. Constructing this surface boils down to solving an unusual boundary value 
problem for the laplacian on the transverse space (fiat in the case of four- dimensional 
collisions). We will follow a similar approach for shocks in anti-de Sitter space. 

Heuristically, a marginally trapped surface S is the limit of a trapped surface where one 
of the null normals propagates inward and the other propagates in a direction that is neither 
inward nor outward, but only forward. To give a more precise and useful definition, consider 
a null basis {i'^,n^) for the normal plane to S at any given point. By convention, i^^ is 
outward pointing and is inward pointing. Both are required to point forward in time. S 
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is defined by the requirement that it is closed and spacehke and that the expansion of £^ 
should vanish: 

^ = /i''-V/, = 0, (97) 

where h^j^y is the induced metric on S. The equation (p7|) is the mathematical expression of 
the heuristic notion that the outward pointing normal propagates neither outward (which 
would lead to positive expansion) nor inward (negative expansion). 

Following [SI [HI Uni [11], we look for a trapped surface S, made up of the union of two 
pieces, S = SiVJ 82- The first piece Si lies in the null hypersurface u = with v < Q while 
the second piece ^2 has v = Q and n < 0. The hypersurfaces Si will be found by looking 
for appropriate codimension 2 surfaces with vanishing expansion ( 1971) . supplemented by the 
boundary condition that the outward pointing normal to S are continuous at the intersection 
C = iSi n S2- A cartoon of this construction is shown in figure H], but it is important that S 
does not exist at a fixed time t: parts of it exist for all times t < ol^ 

Working out the shape of Si and S2 for colliding shock waves in AdS^:) is facilitated by 
a change of coordinates 

V -^v + (t){x\z)Q{u) (98) 

where (f){x\z) is the function appearing in (15^ and ([Ij) and 6(m) is the unit step function 
defined in f l48l) . After this change of coordinates, the components of the metric are everywhere 
finite, but with jump discontinuities at u = 0|lj This means that geodesies have no coordinate 
discontinuities. We define Si as the surface 

u = v = -ipi{x\z) (99) 

for some function ipi defined on the transverse space H£)~2- More precisely. Si is the region of 
the submanifold where ipi > 0, and its boundary C is the curve on H£)_2 (at u = v = 0) 
where ipi = 0. The 0{D — 2) symmetry of a head-on collision means that C must be a 
{D — 3)-sphere, q = qc for some constant qc- An obvious basis for the normal space to Si 
in the cotangent bundle is {du, dv + dipi). So one must be able to express the outward null 
vector as 

i^^Ux^" = Adu + B{dv + dilJi). (100) 



^It is strange indeed to think that entropy exists at times t < 0. The right interpretation is that the 
trapped surface puts a lower bound on the amount of entropy that must eventually be created. 

^"^One can also use a more sophisticated shift, similar in spirit to the one used in |50j where Xi and z are 
shifted and the resultant metric becomes continuous at u = 0. See for example [6]. This leads to the same 
results. 



25 



Since i^f}^ is null (because i^f^^i'^-^''^ = 0), forward (because f > 0) and outward (because 
i"" < and w = is inside the surface), we find 

A = -{dAr B = -^-. (101) 

Here {dipiY niay be equivalently computed using the metric on Hd-2 or the metric on AdS/j. 

By symmetry (assuming that the momenta of the shock waves are equal and opposite), 
^2 must be the image of Si under the interchange of u and v. Thus 

£(.^)dx>' = ~id^lj2)^dv -^{du + dij2) , (102) 

and moreover ipi = 1^2 since the collision is head on, so let us denote them simply as ip. 
Continuity of the outward null normal across C means that £(1) = £(2) when u = V = %jj = 0. 
Thus we require 

(d^Y = ^ on C. (103) 

The Poincare coordinates {u,v,x^,z) do not make the 0{D — 2) symmetry manifest. In 
place of u and v, it is better to use Lu/z and Lv/z, because these combinations are just 
X° — X^~'^ and X° + X^~'^, and the global coordinates X^^ make the 0{D — 2) symmetry 
apparent. Correspondingly we define 

^ = -iIj. (104) 

z 

Evidently, (11031) together with the vanishing of '?/^ on C are equivalent to 

^|c = (9^)2|c = 4. (105) 
Plugging either of fllOOj) or fll02p into WJ\ leads to the condition 

□//._.-^^)(*-<f)=0. (106) 
The most general solution to (11061) respecting the 0(-D — 2) symmetry is 

$ = $ + C$_, (107) 
where C is an integration constant, $ is given as in (l68l) . and $-(g) = 1 + 2g as in (ITHj) . 
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Note that \1/ > inside the curve C because ip is positive. So (11051) becomes 

2L 

vl>(gc)=0 ^>\qc) = -^====. 

^Jqc[l + qc) 

Combining (11071) with (llOSp one obtains an equation for qc- 



(108) 



^\qc 



l + 2qc 

This can be conveniently rewritten as 

W.{qc 



'^{qc 



2L 



y^qcil + qc) 



0. 



(109) 



Jqc{l + qc) 



where 



W4q) = <l>(g)$'„(g) - <l>'(g)<l>_(g) 



(110) 



'111) 



is the Wronskian of $ with Starting from (1771) . it is easily checked that the Wronskian 
satisfies the following first order differential equation: 



q{l + q)d, + ^^{l + 2q) 



= 167rGflL2$„(g)p(g) . 



(112) 



The solution to (ITT^ ifl 
W4q) 

where 



TlG 



D 



[q{l + q)] 



{2-D)/2 



^(g) = 2^-3L^-2(Vo15 



dqo [go(l + go)] 



(D-4)/2 



1 + 2go)p(go) 



(113) 



(114) 



Comparing (I114p to (l83|) . one sees that E{q) can be interpreted as the energy of the massless 
particles inside a radius q in the transverse space i?D-2- Putting (II 10!) together with (I113p . 
we arrive at a simple relation between qc and the energy Eq = E{qc) of massless particles 
inside a radius qc: 



EcGd 



}D-6 



{Vols 



D-3\ 



IT 



;i + 2gc) [qc{l + qc)] 



(D-3)/2 



(115) 



^^There is a more general solution, obtained by adding a multiple of [q{l + q)]^^ ^^^^ to ()113p . But this 
solution is singular at g = 0, and it corresponds to adding an additional, finite-energy massless particle at 
q = 0. 
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47r ^ 



4 v/?c(l + gc)^^ 7r(yr7^-i) 



5 2gc(l + gc)(l + 2gc) 7r(xcy 1 + - sinh"^ xc) 

6 2nqf{l + gc)'/'(l + 2gc) y [(x^ - 2)y'T7^ + 2] 

Ifivr vr^ / 

7 —qUl + gc)'(l + 2gc) -[{2x1 - 3)xcy/i + + Ssinh-^xc] 

Table 1: Values of Ec and ^trapped for 3 < -D < 7, in terms of the chordal distance variable 
qc and Xc = 2^/q^{TTqc)■ 

Table [1] includes values of Eq for 3 < D < 7. 

Once qc is known, it is straightforward to calculate the area of the marginally trapped 
surface S. Because Si is embedded in the null hyperplane u = 0, whose transverse part 
is H£)^2 and whose lightlike direction is parameterized by v, the induced metric on Si is 
identical to the one on H£)_2 inside a radius gcO Thus, the area of Si is just the volume of 
the ball q < qc Hd~2, namely 

= L^^'iVolS^^') r dx^^ . (116) 
Jo y/l + x^ 

In (11161) we have introduced a new radial coordinate 



X 



iL = 2y^^(lTq), (117) 



where r is the radial variable appearing in flM]) . and q is the usual chordal distance variable, 
appearing for example in The area of ^2 is of course the same as of Si. So the entropy 
bound is 

S > 5'ti-appcd = 77^—2^51 . (118) 

The integral (I116p can be expressed in terms of incomplete beta functions, but the explicit 



^^Explicitly: the induced metric on Si is 

£2 



du dv dx^ dx^ dx^ dx^ dz dz 

+ 



di'^diP di°^diP d^'^d^l^ 



where ^" = (^^ ,£,'^,£,^) are any choice of coordinates on Si . The first term drops out because u = identically. 
Choosing ^" = {x^, a;^, z), we immediately recover the metric of H^. 
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form is unenlightening. Table [T] includes values of 5'trapped for 3 < D < 7. 

In principle, qc and xq can be eliminated from the relations (11151) . flll6p . flll7p . and (11181) 
to obtain an explicit dependence of S'trappcd on Eq. In practice, this is difficult to carry out 
explicitly for D > 3. However, when qc ^ l, ii is straightforward to extract the leading 
power law dependence: again for D > 3, 

/ rD-2\ 1/-D-2 

Strapped ^ Co ( j (i^c^) (119) 



with 

Cd 



(22D-7^D-3 Vol 5^-3) ^^20) 

^-3 ■ 

When the shock wave is point-sourced, Ec = E, the total energy. The special D = 5 case of 
(11191) . with Ec = E, was quoted in ([8]). When calculating Strapped in shock wave collisions 
intended for comparison with heavy-ion collisions, it suffices to use the leading power law 
dependence indicated in ( 11191) . The value of S'trapped obtained by using the exact parametric 
relations given in table [1] are only about 0.1% different when qc ~ 38, which is the value 
corresponding to EL ^ 4.3 x 10^. The correction is so small because it is parametrically 
O (^^) for = 50 

The case D = 3 is evidently special. The leading order expansions (11191) and (11201) don't 
work, but instead we can eliminate qc and xc altogether from the expressions in table [1] and 
find 

5'trapped-^ = COsh"^ inGsEc , (121) 

where we must have AirGsEc > 1 to form a trapped surface at all. Let us focus on point- 
sourced shocks, so that Ec = E, the total energy of one shock wave. Remarkably, the 
result (I12ip coincides with a well-known exact result [51] on light-like particles in AdSa. The 
quantities denoted and e in [51] (see for example (2.9) and (3.1) of that paper) should be 
identified as 

pO = tane = 471^3^ . (122) 
Using (4.3) and (4.6) of [51] and noting that in this paper L is set to 1, one finds 



^■^The parametric dependence of the leading correction to (|119p depends on D: it is 0{l/qc) for Z? = 4 
and 0(l/gc) ioT D^6 and £» = 7. 
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which indeed agrees with (11211) . Also, it was shown in [51] that for p° < 1, no horizon 
would form in the collision of the light-like particles; instead, they would merge to form a 
massive particle with no horizon around it. This matches with the observation that one 
needs AttG^E > 1 to form a trapped surface. 

As noted in [51], black hole circumference can't increase except by sudden events like 
its formation, because gravity is non-dynamical in AdSa. So it makes sense that the bound 
S > Strapped should be saturated. A partial translation of this statement to the dual field 
theory is that no entropy increase is possible after formation of a thermal state because there 
is no viscosity. Bulk viscosity is forbidden by conformal invariance, and there is no shear in 
1-1-1 dimensions. Similar observations have recently been made in [15| \5^ . 

So S* = S'trapped in AdSs, at least for point-sourced shocks. As pleasing as this result 
appears, we remain puzzled on a conceptual level: the field theory dual probably enjoys 
some form of integrability and/or holomorphic factorization which permits only forward 
scattering. If one tries to collide pairs of shock waves in the dual field, forward scattering 
says that they should pass right through each other. It would be nice to imagine that after 
the collision, the shocks "bleed" detritus from their leading edge, and this detritus comes to 
rest (on average) as a fully-formed thermal medium. But we could not see how to describe 
such a process hydrodynamically without violating local momentum conservation. It should 
probably be kept in mind that the exact results of [51] rely on the assumption that colliding 
particles mergeo This is certainly a minimal assumption, but it does not seem inescapable. 



5 Discussion 

The scaling S'trapped oc E"^^^ is the most distinctive feature of point-sourced shocks in AdSs 
and, as we have already remarked, it will conflict with data if the S oc E^^"^ behavior observed 
in heavy ion data to date extends significantly above the scale of RHIC collisions, or if the 
slower increase of S with E predicted by CGC calculations is realized. Let's consider how 
we might modify the colliding shocks so as to be consistent with a slower increase of S with 
E. Instead of asking how we could suppress entropy production while keeping energy fixed, 
it is intuitively easier to consider collisions with fixed Strapped and ask how we could add 
energy to them. An obvious approach is to use a halo effect: diffuse energy density can 

-"^"^The additional dynamical claim of [51] is that the holonomy of the end state particle is the product of the 
holonomies of the initial state particles. This product rule explains how the dilference between p° < 1 and 
p° > 1 arises: in the former case, the product of holonomies of initial state particles is a rotation, meaning 
that its fixed set is timclike, corresponding to a massive particle; while in the latter case, the product is a 
boost, corresponding to a spacclike geodesic interpreted as the future singularity inside the black hole. 
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be added outside of where the trapped surface forms, and because it's so diffuse, it doesn't 
cause significant entropy production during the colhsion. But instead of spreading out the 
energy density in R^'^, we're going to spread it out in AdSs in such a way that (T^n) in 



the gauge theory doesn't change. As noted in section I3.3[ there is a large freedom in how 
to do this: any gravitational source in AdSs which has the 0(3) symmetry preserved by a 
massless point particle and which is localized at n = (or v = 0) will give the same (Tmn) 
(but, presumably, with different higher point functions of Tmn and other operators). We will 
consider a particular class of 0(3)-symmetric sources that make it easy to obtain explicit 
formulas using results already established. Namely, let the gravitational stress tensor be of 
the form 0751) where 

for some positive constant a. The overall normalization was chosen so that E is defined as 
in fl83|) . To simplify notation, let us introduce scaled forms of the energy and entropy: 

^ _ 2G^E ^ _ Os'S'trappcd (125) 

It is helpful to keep in mind that E and 5* are numerically large: with G5 and L chosen 
as indicated in section [1], E = 19.7 TeV and S = 35000 translates into E i=s 4.5 x 10^ and 
S ^ 5900. 

From the third entry in table [H for S'trappcd one can show that 

S^Mh (126) 

where the approximate equality becomes more accurate at large S. Using flll4l) -( |TT5l) to- 
gether with (I126p . one arrives at 

^3/2 

E ^ , (127) 

1 - (1 + Vs)- 

where again the approximate equality becomes more accurate at large 5*. Taking a ^ 00 
for fixed S, one recovers from (11271) the leading order result E ^ S^^"^, which is equivalent 
to flll9p for D = 5. This is because in the large a limit, the distribution of energy (11241) 
becomes pointlike in AdSs. 

By allowing a to be a function of S, one can evidently persuade the result (11271) to 
conform to a desired scaling relation, at least for large S: for example, to get E ^ KS"^ for 
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some constant one would choose 



'I" ■ . ' (128) 



\/S\ogS \ 2AV 

at large S. Small a means that the cloud of matter in AdSs is very diffuse, so that most of 
the energy is in the halo where little entropy production occurs. 

The strategy outlined here for obtaining a scaling ^trapped oc E^^"^ seems to us ad hoc. Why 
would we disperse the matter in AdSs as we increase the boost factor? Some guidance from 
other physical principles — perhaps saturation physics — is needed to specify more precisely 
what initial state we should choose in the holographic dual. Another possibility to explain a 
deviation from the E'^^^ scaling is that even slight broadening of the matter distribution in 
the longitudinal direction would lead to a substantial reduction in the production of entropy. 
This possibility is hard to assess without a more careful analysis of trapped surfaces in 
geometries with longitudinal smoothing. We leave such an investigation for future work. Of 
course, we should not entirely neglect the possibility that future heavy-ion data will show a 
markedly faster increase in total multiplicity with energy than the Landau model predicts. 
If flTB]) turns out to be about right, we would see it as evidence that the trapped surface 
computation captures an important aspect of the overall dynamics of the collision. 

Although entropy estimates have been our main focus, the 0(3) symmetry that we noted 
in head-on collisions of point-sourced shocks has independent interest. In this context, it is 
interesting to note the proposal of [S3], according to which the saturation scale Qg should 
vary across the transverse plane in a fashion similar to 



Q'(-\^) = JT^J^f^Q''r, (129) 

where L is taken to be the transverse size of the hadron under consideration!^ The form 
(11291) was proposed in order to economically accommodate the known power law behavior 
Qs ~ l/[(x^)^ + (x^)^] at large transverse x together with a finite maximum at = = 0. 
It is assumed to arise from feeding an initial state (the hadron at rest) which explicitly 
breaks conformal symmetry into evolution equations (BFKL and generalizations) which are 



^^The saturation scale is the typical transverse momentum of color field configurations in a highly boosted 
nucleus. It arises when perturbative splitting causes the phase space density of gluons to become of order 
l/tts, at which point recombination of gluons cannot be neglected. For gold-gold collisions at RHIC, an 
approximate value for Qs is 1 GcV. Sec for example }54j for an introductory account of the saturation scale 
and related ideas. 
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conformally invariant, at least in the leading-log approximation. 

The form (11291) is tantalizingly similar to fl86p . hinting that the initial state and early 
dynamics of the collision might be closer to respecting the 0(3) symmetry we have found 
than one would expect from comparing the shapes of the Woods-Saxon and AdS profiles 
shown in figure HI Indeed, the profile fl86|) can be recognized as special even without recourse 
to the gauge-string duality: it is uniquely specified, up to the choice of E and L, by invariance 
under the 0(3) remnant of the conformal group that we mentioned in section [1] and identified 
explicitly at the end of section 13.11 We speculate that this 0(3) may be approximately 
realized in central heavy ion collisions. If it were realized exactly in the initial state, it would 
be preserved during the collision to the extent that the relevant dynamics — perturbative or 
strongly coupled — respects conformal symmetry. 

The 0(3) symmetry has a particularly simple realization on the S*^ x R boundary of global 
AdSs. To make it explicit, let's introduce explicit polar coordinates in AdSs as follows: 

= Vp^Tl2 cos ]_ x° = + ^2 sin ]_ 

Ij Lj 

= psin'?/;sin'!9cos0 = p sin sin ■i? sin (130) 

= p cos i^) X^ = p sin i\) cos . 

The angle is the usual azimuthal angle around the beam-line, but the other angles do not 
have such a familiar interpretation: in particular, is not the angle relative to the beam. 
The boundary metric is 

ds^ = -dr^ + L'^idip^ + sin^ d^^ + sin^ ^jj sin^ Mcp'^) . (131) 



As we saw in the discussion at the end of section [3TT| it takes a massless particle a global time 
At = ttL to traverse AdSs. Let's say that a right-moving particle starts on the boundary 
at = 7r at global time r = —ttL/2, and a left-moving particle starts at the same time at 
ip = 0. The propagation of these particles toward one another is dual to an expansion of 
light-like, spherically symmetric shock waves from the insertion points aX ip = tt and 0. At 
times —TiL/2 < r < 0, the stress tensor on S*^ x R is 



1 ~ E 
(^-) = ;^(^'AV^) = 4^^3sin2^ 

E 

AttL^ sm ip 



(132) 
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with other components vanishing. The first term in square brackets of each of the exphcit 
expressions in (11321) is due to the right-moving shock, and the second term to the left-moving 
shock. The notation Tab for the stress tensor on 5*^ x R reminds us not only that there is 
a non-trivial conformal mapping between this form and fl86l) . but also that the stress tensor 
picks up an anomalous vacuum contribution, proportional to the metric, in the course of this 
mapping, which we have excluded from Tab- 
It would be interesting to start from (11321) . or some alteration of it where the delta 
functions are softened into sharply peaked but smooth functions, and evolve it forward with 
hydrodynamics, or with some combination of a heuristic treatment of thermalization for 
< r < Ttherm foUowcd by hydrodynamics. Although such evolution would be a purely 
one-dimensional problem, it would incorporate a combination of radial and longitudinal flow 
in the original Minkowski-space conformal frame. 

It is worth noting that an 0(2) remnant of the conformal group is preserved even when 
the collision is not head-on. (If D > 5, then this would be an 0{D — 3) symmetry.) This 
is easiest for us to see by considering test particles in AdSs. In Poincare coordinates, the 
right-moving particle travels on a trajectory with 



X 



3-t = 3.2^0 z = L, (133) 



and the left-moving particle's trajectory is 



X 



3 



-t x^=0 x^ = b z = L, (134) 



where b is the impact parameter. In global coordinates, these trajectories take the form 

right-moving: = X° X^ = X"^ = X^ = X'^ = L 

{X^ = -X° = X^ = b (135) 

X^ = -— X-' = L + — . 
2L 2L 

The total configuration (11351) is invariant under the full 0(3) subgroup iff & = 0, but for 
6 7^ it is still invariant under the 0(2) subgroup that preserves the vector {X^, X"^, X*) = 
(0,6, —b^/2L). Because the dynamics respects the full conformal symmetry, the final state 
should be invariant under 0(3) or 0(2), according as 6 = or 6 7^ 0. 

Clearly, estimating total entropy production is only one facet of describing colliding 
shocks in anti-de Sitter space. Ideally, one would like to understand the process of thermal- 
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ization and the subsequent hydrodynamical flow. In the case of central colhsions, imposing 
the 0(3) symmetry means that the gravity calculations are effectively 2 + 1-dimensional. 
An optimistic view is that a fairly full account could be made by somehow matching a 
marginally trapped surface computation to a late-time description that fuses a linearized 
treatment of non- hydrodynamical quasi-normal modes, as described in [551 [56], with a non- 
linear treatment of hydrodynamical modes, as described in ^7\. In such an account, the 
local thermalization time might be expected to be several times the relaxation time of the 
non-hydrodynamical modes: in total, roughly 0.3fm/c for conditions comparable to central 
RHIC collisions, according to the estimate of [56]. An alternative approach to describing 
thermalization following a collision of two shocks in AdSs has been suggested in [B]. 
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